A New Method to Search for CP violation in the Neutrino Sector 
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New low-cost, high-power proton cyclotrons open the opportunity for a novel precision search for 
CP violation in the neutrino sector. The accelerators can produce decay-at-rest neutrino beams 
located at multiple distances from a Gd-doped ultra-large water Cerenkov detector in order to 
search for CP violation in P M — > at short baseline. This new type of search complements presently 
proposed experiments, providing measurements that could lead to a substantially better exploration 
of CP violation in the neutrino sector. 
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I. INTRODUCTION 

With the discovery of neutrino oscillations, particle 
physicists have been inspired to develop theories that 
explain very light neutrino masses. The most popular 
models invokes GUT-scale Majorana partners which can 
decay, producing a matter-antimatter asymmetry in the 
early universe through the mechanism of CP violation. 
Observation of CP violation in the light neutrino sector 
would be a strong hint that this theory is correct. 

To incorporate CP violation, the light- neutrino mix- 
ing matrix is expanded to include a CP violating phase, 
5c p- Sensitivity to Sep comes through muon-to-electron 
flavor oscillations at the atmospheric mass-squared dif- 
ference, Am^. The oscillation probability, neglecting 
matter effects, is given by pQ: 

P = sin 2 6> 23 sin 2 26» 13 sin 2 A i3 
=F sin 5 cp sin 2f?i 3 sin 2# 23 sin 26*i 2 

x sin 2 A13 sin A12 
+ cos S cp sin 26>i3 sin 20 23 sin 2#i 2 

x sin A 13 cos A 13 sin A 12 
+ cos 2 6> 23 sin 2 20 12 sin 2 A 12 , (1) 

where A^ = AmfjL/AE^, and — (+) refers to neutrinos 
(antineutrinos) . 

In Eq. [T] aside from 8c p, all but two of the parame- 
ters (#i3 and the sign of Arn^j) are well known. Further 
precision on these known parameters is expected in the 
near future (see Table [I]). With respect to #13, global fits 
report a non-zero value at the ~ la level This pa- 

rameter drives the amplitude for the CP violating terms 
in Eq. [I] and therefore sets the level of technical difficulty 
for observing CP violation. The unknown sign of Am^, 
referred to as "the mass hierarchy," affects the sign of 
term 3 in Eq. [T] Matter effects modify Eq. [T] and are 
sensitive to the mass hierarchy. 

We propose a new method to search for CP violation 
by comparing absolute neutrino rates in a single detector 
that is illuminated by neutrino beam sources at multiple 
distances. This method explores the Sep dependence 
inherent in Eq. [I] in a different way and with comparable 
precision than the planned long-baseline program. 



The neutrino sources would be based on commercially- 
developed, small (2.5 m diameter), high-power proton cy- 
clotron accelerators that are under development [7|. A 
250 MeV, 1 mA cyclotron is under construction at MIT 
and a GeV-energy, megawatt-class cyclotron is presently 
under design. When in production, because of new, inex- 
pensive superconducting technology, these machines are 
expected to cost 5% of a conventional proton accelerator 
(< $20M). These machines will demonstrate the beam 
physics and engineering needed for this experiment. 

This experiment requires accelerators that target 2 
GeV protons at 2.5 mA during a 100 /is pulse every 
500 /is, delivering 9.4 x I0 22 protons per year to a beam 
stop. The result is a high-intensity, isotropic, decay-at- 
rest (DAR) neutrino beam arising from the stopped pion 
decay chain: 7r + — > + /i + followed by /i + — > e^v^v^. 
The flux, shown in Fig. [t] has an endpoint of 52.8 MeV. 
Each cyclotron produces 4 x 10 22 /flavor/year of 24, 
and P„. The v e fraction in the beam is very low (~ I0~ 4 ) 
because most ir~ are captured before decay. 

This experiment utilizes an ultra-large water Cerenkov 
detector such as Hyper-K [5], MEMPHYS [9], or the de- 
tector proposed for the Deep Underground Science and 
Engineering Laboratory (DUSEL) [10]. We use a 300 
kton DUSEL detector as our model. Water provides a 
target of free protons for the inverse beta decay (IBD) 
interaction: v e + p — > n + e + . 

IBD interactions arc identified via a coincidence signal. 
The first signal is from the Cerenkov ring produced by 
the positron. The second signal is from capture of the 
neutron. The signal from neutron capture on protons, 
which produces only a single 2.2 MeV 7, is too feeble to 
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FIG. 1: Energy distribution of neutrinos in a DAR beam 
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FIG. 2: Proposed experimental layout. Each phase corre- 
sponds to a 5-year data-taking period with each of the three 
source sites on for 20% of the time. (Phase 2 is an example; 
the number of cyclotrons will be optimized based on Phase 1 
measurements.) 



be efficiently observed in a large Cerenkov detector. For 
that reason, doping with gadolinium (Gd), which has a 
high capture rate and a short capture time, is proposed 
[llj . Neutron capture on Gd produces multiple photons 
totaling ~8 MeV which can be observed with 67% ef- 
ficiency |12j . Gd-doping of water is under development 
13J . In IBD events, the energy of the neutrino is related 
to the energy of the positron and kinetic energy of the 
neutron, K, by Ev = E e + + (M n - M p ) + K [14], 

The accelerators are placed at three separate baselines, 
L, for the — > v e search (see Fig. [2|. Each set of ac- 
celerators at the three sites runs for a 20% duty factor 
interspersed in time; the other 40% of the time is used to 
collect beam-off data. The longest baseline is L = 20 km, 
where ~40 MeV are at oscillation maximum for the 
measured value of Ato^. In this case, term 2 of Eq. [l] 
will contribute while term 3 will not. A mid-baseline ac- 
celerator site is located at L = 8 km, where term 3 of 
Eq. [T] is non-zero. The shortest baseline accelerator site 
is at L =1.5 km, directly above the DUSEL detector. 
At this location, the neutrino-electron elastic scattering 
rates are sufficient to allow precision measurement of the 
neutrino flux normalization as discussed below. We pro- 
pose multiple accelerators at the various locations (see 
Fig. [2J , which is feasible because of the low cost per ma- 
chine. For the L of a given event to be well-determined, 
the beams from each location are staggered in time. 

Term 3 of Eq. 1 leads to an inherent ambiguity be- 
tween Sep and the mass hierarchy. We assume that on 
the timescale of the measurements presented here, the 
mass hierarchy will be measured using LBNEpjTj or at- 
mospheric neutrino measurements. [T5] . For the sensitiv- 
ity estimates, we choose the normal hierarchy as the ex- 
ample model. 

We present a two-phase experiment, where the design 
has flexibility in the second phase. In our case, each 
phase represents five years. The first phase explores the 
oscillation space at ^3cr and requires one accelerator at 
the near location, two at the mid-location and three at 
the far location. Once a signal for CP violation has been 
localized, the strategy for the second phase can be de- 
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Amji x 10~ b eV 2 


7.65 


0.23 




7.65 






Am§! x lCT 3 eV 2 


2.40 


0.12 




2.40 


0.02 
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sin 2 (26»i2 ) 
sin 2 (2<9 23 ) 
sin 2 (2f9 13 ) 


0.846 


0.033 




0.846 






1.00 


0.02 
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1.00 


0.005 


A 


0.06 


0.04 


m 


0.05 


0.005 
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TABLE I: Left: Present values and uncertainties for oscilla- 
tion parameters, reported in the references. Right: Future 
expectations used in this study. 



termined. As an example for phase 2, we have chosen a 
design with one, two and seven accelerators. This design 
matches the sensitivity of a run of LBNE with 30 x 10 20 
protons on target in neutrino mode followed by 30 x 10 20 
protons on target in antineutrino mode |10) . 

This design is unique among proposals for CP viola- 
tion searches and complementary to the present plans 
[2]] . The measurement is done with antineutrinos, 
while all existing proposals rely heavily on neutrino data. 
The neutrino production systematics are different and 
well-controlled. Because of the low energy, the inter- 
action systematics are also different. Varying L, while 
employing a single detector is novel and reduces system- 
atics. A two-phase program which allows an optimized 
measurement strategy is powerful and potentially cost- 
saving. 



II. EVENT TYPES AND BACKGROUNDS 

The energy range of the analysis is 20 < E v < 55 
MeV. The lower cut renders potential backgrounds from 
radioactive decay and spallation negligible. The upper 
cut is chosen because of the 52.8 MeV signal endpoint. 

In this energy range, three types of interactions must 
be considered. First is the IBD signal, with an estimated 
reconstruction efficiency of e re con — 67%, based on stud- 
ies for Super-K [12 . Second is v e + O ->• e" +F (V e O"), 
with a suppressed cross section relative to IBD due to 
nuclear effects [161 117) and no neutron capture. Third 
is neutrino-electron elastic scattering (z/ e ES) v e + e~ — ¥ 
v e + e~. This is separable from v e O by angular cuts 
[T51 Hi?] and, so, can be used for the flux normalization. 
For v e O and f e ES events, we estimate e recon = 75% [18] . 

Beam-off backgrounds arise from three sources: atmo- 
spheric v^p scatters with muons below Cerenkov thresh- 
old which stop and decay ("invisible muons"), atmo- 
spheric IBD events, and diffuse supernova neutrinos. 
These are all examples of correlated backgrounds; acci- 
dental beam-off backgrounds are estimated to be negligi- 
ble. The rates of these correlated backgrounds are scaled 
from analyses for the GADZOOKS experiment [TT]. The 
interaction rates of the beam-off backgrounds are well- 
measured during the 40% beam-off running fraction. 

The beam-on backgrounds have both accidental and 
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Event Type 


1.5 km 8 km 20 km 


IBD Oscillation Events (E„ > 20 MeV) 








8c p = 0°, Normal Hierarchy 


763 


1270 


1215 


" , Inverted Hierarchy 


452 


820 


1179 


8c p = 90°, Normal Hierarchy 


628 


1220 


1625 


" , Inverted Hierarchy 


628 


1220 


1642 


S C p = 180°, Normal Hierarchy 


452 


818 


1169 


" , Inverted Hierarchy 


764 


1272 


1225 


5 C p = 270°, Normal Hierarchy 


588 


870 


756 


" , Inverted Hierarchy 


588 


870 


766 


IBD from Intrinsic V e (E„ > 20 MeV) 


600 


42 


17 


IBD Non-Beam (E„ > 20 MeV) 








atmospheric v^p "invisible muons" 


270 


270 


270 


atmospheric IBD 


55 


55 


55 


diffuse SN neutrinos 


23 


23 


23 


u e -e Elastic (E„ > 10 MeV) 


16750 


1178 


470 


v e - Oxygen (E„ > 20 MeV) 


101218 


7116 


2840 



TABLE II: Event samples for the combined two-phase run for 
sin 2 20i3 = 0.05 and parameters from Table [I] (future). 



correlated sources. The accidental backgrounds arise 
from the v e in the beam which are followed by a neutron- 
like event. This background is estimated to be very small 
using the measurements from the Super-K Gd-doping 
study[12]. Correlated backgrounds are produced by the 
intrinsic v e content of the beam. This is reduced by care- 
ful design of the beam stop. We assume a long water tar- 
get embedded in a copper absorber, surrounded by steel 
for each cyclotron. We target at 100 degrees from the de- 
tector to reduce decay-in-flight (DIF) backgrounds. For 
2 GeV protons on target, the pion yields, n + /p (n~ /p), 
are 0.43 (0.056). with a tt+ (it-) DIF fraction of 0.017 
(0.025). The decays of ^ _ 's are further reduced by a fac- 
tor of 8.3 since the [i~ will be captured before decaying. 
The compact design decreases the D e production com- 
pared to LAMPF/LANCE by a factor of 1.7. The beam 
has a v e /v e ratio of 4 x 10~ 4 . 
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FIG. 3: Phase 1 and Combined-phase sensitivity to #13 7^ 
at 3cr and 5a. 
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FIG. 4: Correlated sensitivity (1 and 2a contours) to Sep and 
sin 2 2#i3 for the combined two-phase running. 



III. SYSTEMATIC ERRORS 

The multiple-DAR-source design is an elegant choice 
for precision oscillation physics because the beam and 
detector systematics are low. The shape of the DAR 
flux with energy is known to high precision and is com- 
mon among the various distances, thus shape compar- 
isons will have small uncertainties. The neutrino flux 
from the three distances is accurately determined from 
the direct measurement of the 7r + production rate using 
^ e ES events from the near accelerator. Existing proton 
rate monitors assure that relative proton intensities are 
understood to 0.1%. The interaction and detector sys- 
tematic errors are low since all events are detected in a 
single detector. The IBD cross section for the signal is 
well-known [14]. The fiducial volume error on the IBD 
events is also small due to the extreme volume-to-surface- 
area ratio of the ultra-large detector. 



The largest systematic errors arise from the i/ e ES sam- 
ple used to determine the absolute normalization of the 
flux. The error on this cross section is 0.5% due to a 0.7% 
uncertainty from the NuTeV sin 2 9w measurement [20] . 
We assume a 2.1% energy scale error [5T] which leads to 
a 1% error on the DAR flux when a E vis > 10 MeV cut 
on the events is applied. The v e events on oxygen and 
IBD events with a missing neutron can be separated from 
the f e ES sample since the angular distribution of f e ES 
events is very forward-peaked, while these backgrounds 
have a broad distribution [TJJ]. We take the uncertainty 
from contamination by these events to be negligible. The 
uncertainty of the electron-to-free-proton ratio in water is 
also very small. Adding the ^ e ES systematics in quadra- 
ture, the systematic error on the IBD flux is expected to 
be 1.1%. For the total error on the IBD flux, one must 
then add the z^ e ES statistical error in quadrature, which 
depends on the running period. 
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sin" 26»i3 /Scp 


-180 


-90 





90 135 


0.01 


52.5 


47.2 


29.0 


38.8 48.5 


0.05 


21.6 


21.9 


19.5 


25.6 30.4 


0.09 


18.3 


19.9 


17.6 


23.8 26.3 


0.01 


51.3 


45.5 


26.9 


36.7 46.9 


0.05 


19.9 


21.2 


18.2 


24.2 29.6 


0.09 


16.8 


19.2 


16.5 


22.6 25.3 



TABLE III: The la measurement uncertainty on Scp for var- 
ious values of sin 2 #13 for the combined two-phase data. Top: 
Systematic and statistical errors. Bottom: Statistics only. 



that of LBNE, but is inverted in its 8c p dependence |10j . 

The combined two-phase running yields a 4.1o mea- 
surement of Sep at the test poi nt o f sin 2 2#i3 = 0.05 and 
90°, as shown in Table 



5c 



III 



(top). The correlated 



measurement uncertainties are shown in Fig. [4] for 1 and 
2cr contours. Table III (bottom), which provides the sta- 
tistical uncertainty, indicates that the measurement is 
statistics- limited. 



The other significant systematic error is on the effi- 
ciency for neutron detection in IBD events. To reduce 
uncertainties, neutrons are tagged via timing rather than 
position reconstruction. This leads to an inefficiency for 
neutrons outside of the timing window with a systematic 
uncertainty of 0.5% [12] . 



IV. OSCILLATION SENSITIVITIES 

Sensitivity estimates were made by calculating the x 2 
for a given set of predicted events and investigating the 
X minimization and excursion, using a method simi- 
lar to Ref. [5]. Matter effects were included in the fit, 
but are negligible due to the short baseline. System- 
atic uncertainties (see Sec. Ill and Table IT| were con- 

of rom/an rvi r n I I i ortn r>An F i'i hnti i~m o T"*-\ i It ^ -\ at t \~t TArm 



strained by pull-term contributions to the x of the form 
(ki — l) 2 /of, where of are the uncertainties. Results are 
given for two scenarios (see Fig. [2| : Phase 1 and Phase 
1+2 combined. Table |H| gives the event samples associ- 
ated with the combined phases, for the various classes of 
oscillation, background, and calibration events. 

The sensitivity for observing a non-zero value for #13 at 
the 3 and 5a CL as a function of Sep is shown in Fig. [3] 
for combined two-phase running. This sensitivity meets 



V. CONCLUSIONS 

We have described a novel experiment to search for 
CP violation in the neutrino sector. This experiment 
is relatively low cost, employing multiple high-powered 
compact cyclotrons located at L — 1.5 km, 8 km and 
20 km from a large water Cerenkov detector. Using the 
example of the DUSEL detector, for sin 2 29^ = 0.05, the 
CP violation parameter Sep can be measured to the level 
of LBNE (> 4o-) in a 2-phase 10- year run. 

The complementary nature of this measurement makes 
it a compelling addition to the program. This experiment 
will probe for CP violation with antineutrinos, unlike 
the present program which relies heavily on neutrino in- 
teractions. The systematics are also quite different and 
low compared to present planned experiments. As a re- 
sult, this experiment will provide a powerful input to our 
global search for new physics in the neutrino sector. 
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